Hypothesis 1: Thermosaline niche shifts

The first hypothesis formulated for this chapter predicts pike to differ significantly in thermosaline niche between life stages, with young pike inhabiting a warmer, less saline niche that continuously becomes colder and more saline as pike grow. To analyze this, I calculated the mean d18O and Sr/Ca values for each individual and year. I assigned an area based on whether the fish was caught in the Bodden or in freshwater. I did not separate the anadromous fish for this step, as there are too few of them to form confident means. The d18O and Sr signatures grouped by years provide the hypothesis with evidence, as the d18O signature of the first year seems to be significantly lower than all consecutive years for marine pike (Fig. 1 A), and the d18O signature increases during the entire life, indicating a constant shift towards colder water during the life history. This also shows for freswater-caught pike (including the anadromous sample, Fig. 1 A), however, the d18O signature for these is significantly lower than for Boddenpike, likely caused by higher d18O water in the lagoons. The Sr/Ca data of the pike also provides evidence for the hypothesis, as the lagoon pike increase in Sr/Ca signal constantly during their life (Fig. 1 B). However, the mean Sr/Ca values seem to reach a plateau after year 5 (if this corresponds to physiological limits of a freshwater fish cannot be answered with confidence in this analysis) and no significant increase occurs thereafter. As expected, the Sr/Ca values of the freshwater pike show no significant increase over the life history. The slight trend that can be observed in Fig. 1 B is due to the anadromous fish included in the freshwater sample, which also increase in salinity in later life stages.

Figure1: Lifelong thermosaline niche of all pike in the sample. A: Mean d18O values for each year of life, B: Mean Sr values for each year of life.Colors indicate life history strategy, marine Bodden pike (blue), anadroous pike (red) and freshwater pike (green)

Figure1: Lifelong thermosaline niche of all pike in the sample. A: Mean d18O values for each year of life, B: Mean Sr values for each year of life.Colors indicate life history strategy, marine Bodden pike (blue), anadroous pike (red) and freshwater pike (green)

Hypothesis 2: Effects of salinity on growth

The first hypothesis “Salinities exceeding a reported optimum of 2 PSU have a negative effect on growth performance of pike that reverses with increasing age and body size of the fish” can only be answered by calculating absolute salinity from otolith Sr/Ca or d18O (note that this “optimum” stands on a rather thin evidence base - one old paper - and thus, the hypothesis could be reformulated to be more general, e.g. salinity bad for small but good for large fish). However, Fig. 2 shows what could be the answer to that hypothesis. The incement widths of all pike from the Bodden, including the anadromous pike, were plotted against the raw Sr/Ca ratio. Pike increment widths seemed to have a negative correlation with Sr/Ca for the first three years as well as for the fifth year, after that, increment widths correlated positively with Sr/Ca. This would agree with the hypothesis and ecological theory, as salinity is less of a stressor the larger a pike becomes, and other factors that might be beneficial for growth, such as marine prey, might outweigh the osmoregulatory investment in later years. As you can see, the \(R^{2}\) values indicate a bad fit for the regression lines. This is likely a result of the differences in SIMS-transects positions and -length, along which the aging has been performed. The transect lines were not chosen in a standardized manner, but rather to maximise the data quality. Therefore, the increment widths are only comparable within an individual. A possible solution for his could be to reassign annuli along a standardized line and assign the yearly mean d18O-values calculated along the transect lines to the new annuli. Another solution could be to calculate an individual metric for growth performance that takes this into account, e.g. \(\Omega\), and correlate that with Sr OR absolute salinity. A third option would be to back-calculate the total length of the fish via a body-proportional equation (e.g. Fraser-Lee), as this would take into account the relative proportions of an individual otolith (ind. radius divided by ind. otolith radius).

Figure 2: Effects of salinity as approximated by mean Sr value per year on growth, as approximated by otolith increment width, in different life stages (years) of pike. Red points are individuals, blue line is a linear regression line, shaded area are the confidence intervals of the regression. Equation and R-square value are annotated in the top left corner of each plot

Figure 2: Effects of salinity as approximated by mean Sr value per year on growth, as approximated by otolith increment width, in different life stages (years) of pike. Red points are individuals, blue line is a linear regression line, shaded area are the confidence intervals of the regression. Equation and R-square value are annotated in the top left corner of each plot

Hypothesis 3

Although not included in earlier plans of this work, it seems a logical step to also include a specific hypothesis for temperature alongside the salinity hypothesis. From the literature, one would expect for temperature to have a positive effect on the growth of young and small fish, but a negative effect on growth of older and larger fish. Plotting the increment widths against mean d18O values, which were assumed a good proxy for temperature, for each year and individual, revealed a negative correlation between d18O and increment width for all years except te oldest fish analyzed (year 10, Fig. 3). As more negative d18O values indicate a lighter signal and therefore higher temperature, this suggests a positive effect of temperature on pike growth throughout nearly the entire age gradient (Fig. 3). However, the fit of the lines, again, is rather bad, for reasons and possible remedies, see above…

Figure 3: Effects of temperature, as approximated by mean d18O value on growth, as approximated by otolith increment width, in different life stages (years) of pike. Red points are individuals, blue line is a linear regression line, shaded area are the confidence intervals of the regression. Equation and R-square value are annotated in the top left corner of each plot. More negative values of d18O indicate warmer temperatures, more positive values indicate colder temperatures

Figure 3: Effects of temperature, as approximated by mean d18O value on growth, as approximated by otolith increment width, in different life stages (years) of pike. Red points are individuals, blue line is a linear regression line, shaded area are the confidence intervals of the regression. Equation and R-square value are annotated in the top left corner of each plot. More negative values of d18O indicate warmer temperatures, more positive values indicate colder temperatures

These exploratory steps hint at potential effects that could be worthwhile to analyze further. Both hypotheses can only be answered conclusively by converting the trace elemental data and isotope values into absolute salinity and temperature units in my opinion. Figs. 2 and 3 might already show an outline of a possible way to analyze this, by extracting the slope from each year class increment vs salinity/temperature model and analyzing the slopes with a GAM (Finger et al., 2016) or random slope approach (e.g. Dingemanse et al., 2010), to determine if the effects of thermosaline niche on growth change significantly throughout pike life history. Another promising approach might be to analyze individual growth patterns, extracting growth performance per year for each individual and modelling growth performance with experienced thermal/saline niche using mixed linear models.

Ecotypes?

Based on the notion of (visually) different plots (as shown exemplary in figure X), I attempted to discern 4 different ecotypes/life strategies from the d18O & Sr data. I aimed at describing a freshwater resident (F-type, Fig. 4 A), an anadromous ecotype (A-type, defined as originating from freshwater but migrating into brackish salinities in later life, Fig. 4 B), a high salinity ecotype (S-type, defined from a high Sr/CA value in the first year and above average Sr-values for the rest of the life, Fig. 4 D), and an intermediary ecotype (I-type, defined from a lower Sr/Ca-value in the first years that´s higher than the freshwater value but lower than Bodden average, and a constant increase in Sr/Ca throughout life, Fig. 4 C). Additionally, based on the d18O value in the first year and the lifelong mean d18O, I assigned three temperature types, warm (W), low d18O in first year and lifelong below global average, intermediary(I), close to global mean in first year and lifelong, and cold (C), above global mean in fist year and lifelong.
Figure 4: Visual ecotypes for pike in the study area. The left plots shows Sr/Ca values during the entire transect from core (birth) to outer margin (point of capture), indicating saline history. The right plots show the d18O values from core to margin, indicating thermal history. A: Freshwater pike, low Sr during the entire life, d18O signature well below global mean during the first year and lifelong, indicating a warm thermal type. B: anadromous pike migrating regularly into freshwater and back into the Bodden. The d18O signature of this pike is below the global mean during the entire life and especially during the first year cycle, indicating a warm thermal type. C: Intermediary saline type, characterized by constant shift towards more saline water during the life history. The d18O signatue is near the gloal mean during the entire life history, indocating an intermediary thermal type. D: High saline ecotype with origin in high Sr and high signal during the entire life, indicating a high salinity ecotype. The d18O signature is above the global mean during the entire life, indicating a cold thermal type.

Figure 4: Visual ecotypes for pike in the study area. The left plots shows Sr/Ca values during the entire transect from core (birth) to outer margin (point of capture), indicating saline history. The right plots show the d18O values from core to margin, indicating thermal history. A: Freshwater pike, low Sr during the entire life, d18O signature well below global mean during the first year and lifelong, indicating a warm thermal type. B: anadromous pike migrating regularly into freshwater and back into the Bodden. The d18O signature of this pike is below the global mean during the entire life and especially during the first year cycle, indicating a warm thermal type. C: Intermediary saline type, characterized by constant shift towards more saline water during the life history. The d18O signatue is near the gloal mean during the entire life history, indocating an intermediary thermal type. D: High saline ecotype with origin in high Sr and high signal during the entire life, indicating a high salinity ecotype. The d18O signature is above the global mean during the entire life, indicating a cold thermal type.

Plotting the measured increments on the otoliths for the different ecotypes shows differences in growth performance in Fig. 5, but they do not appear to be significant within most year classes. Significant differences in later years might be caused by low sample size. The freshwater resident fish show relatively high increment widths, a proxy for growth performance in a given year class, throughout their lifes, except for the first year, however, this ecotype class only consists of 7 individuals. A high relative growth performance was found in the intermediary type in the first year, with medium growth for the rest of its life history. The anadromous ecotype shows high variability in increment widths in the first few years, implying effects other than the general life strategy (such as effects of the stream it was born in or the Bodden it grew up in), however, it was apparent that this type shows a high growth performance in the first three years of life and a lower growth performance in later years. The high salinity ecotype shows a relatively low growth performance in the first three years, a medium growth in year four to six and a relatively high growth in the later years. While there seem to be differences in growth, none of them are significant within a year (except for ages older than 7 years, but note that the sample size become small for old fish). Looking at lifelong growth performance at the individual level might reveal more of these trends and allow to separate spatial effects of Bodden and origin tributary from effects of life strategy/ecotype.

Figure 5: Increment widths of pike assigned to the four ecotypes during their life history

Figure 5: Increment widths of pike assigned to the four ecotypes during their life history

Fig. 6 shows the mean increment widths of pike assigned to three temperature classes, with mean increment width in um on the y-axis and age in years on the x-axis. Addtionally to the ecotypes, which are primarily motivated by the Sr/Ca data, the temperature classes were assigned based on the d18O-signature in the first year and the mean lifelong d18O relative to the global mean of all pike. The three classes cold (C), intermediate (I), and warm (W) also showed differences in growth. The cold temperature type showed lower mean increment widths up until the fifth year of life, with the warm type scoring higher growth particularly in the first three years. In later years, the cold and intermediary type show a higher growth than the warm type. However, the d18O data were not cleaned from the effect of d18O water, so effects of different salinity depending on Bodden/tributary of origin can not be excluded.

Figure 6: Increment width of pike assigned to the three thermal types duing their life history.

Figure 6: Increment width of pike assigned to the three thermal types duing their life history.

Ecological theory would predict a pike to seek out a cooler thermal niche when occupying a high saline niche. The qualitative ecotypes/temperature classes support that notion. Fig. 7 shows the frequency of qualitative temperature classes per saline ecotype. As one would expect, the anadromous and freshwater ecotype mostly consist of the warm and intermediate temperature types, with a small percentage of cold temperature types in the anadromous ecotype. The intermediary saline ecotype consists of mostly intermediate temperature types. The saline ecotype consists of predominantly cold and intermediary temperature types.

Figure 7: Frequency of temperature types among salinity ecotypes

Figure 7: Frequency of temperature types among salinity ecotypes

Looking at clustering

So, in the manner of all exploratory analyses with ecology in mind, of course I had to calculate a PCA with the data. The group colouring in the two plots below are the assigned ecotypes, but I´ve tried many. As you can see, in the first year (Fig. 8), most of the pike differentiate along the d18O and Sr Eigenvectors, with some variance along Na (which is also a proxy for salinity if I remember correctly) and Mg (which strongly correlates with Sr in aragonite). In later years (Fig. 9, I´ve chosen year 5, as this seems to be the year after which the effect of salinity reverses, see above), there´s quite some variance in the data, and the intermediary, high saline and anadromous fish form one big “cloud”, showing that the visual ecotypes might not be the best way to characterize them. However, I´ve also tried to colour them with their respective waterbody of origin in Fig. 10 (or, more precisely, the waterbody they´ve been captured in, we know how far those fish can and will move sometimes). There might be some potential for geolocation of fish based on elemental fingerprints in this data, but since this will not be the main focus of my second chapter, I´ve kept this analysis short and didn´t dive too deep into it just yet.

Figure 8: PCA-biplot showing the PCA for first year pike ordered according to the first two principal components. Colored arrows show the Eigenvectors of the respective elemental data, different colors in the data points show visually assigned ecotypes

Figure 8: PCA-biplot showing the PCA for first year pike ordered according to the first two principal components. Colored arrows show the Eigenvectors of the respective elemental data, different colors in the data points show visually assigned ecotypes

Figure 9: PCA-biplot showing the PCA for fifth year pike ordered according to the first two principal components. Colored arrows show the Eigenvectors of the respective elemental data, different colors in the data points show visually assigned ecotypes

Figure 9: PCA-biplot showing the PCA for fifth year pike ordered according to the first two principal components. Colored arrows show the Eigenvectors of the respective elemental data, different colors in the data points show visually assigned ecotypes

Figure 10: PCA-biplot showing the PCA for fifth year pike ordered according to the first two principal components. Colored arrows show the Eigenvectors of the respective elemental data, different colors in the data points show waterbody of capture of individual pike

Figure 10: PCA-biplot showing the PCA for fifth year pike ordered according to the first two principal components. Colored arrows show the Eigenvectors of the respective elemental data, different colors in the data points show waterbody of capture of individual pike

I also attempted to form discrete clusters. From Fig. 11 below you can see that having four clusters seems to be a good compromise between reducing within group sum of squares and keeping complexity low. These four clusters are shown in Fig. 12. Assigning the clusters to the Sr and d18O Biplot in Fig. 13 yields some insights. There seem to be two distinct freshwater cluster (clusters 3 and 4), whereas the two clusters 1 and 2 seem to be predominantly Bodden clusters. The anadromous fish are all over the place and present in all clusters, showing again that this characterization could use some tweaking. If you have any ideas on how to improve on this, please let me know.
Figure 11: Plot of within groups sum of squares depending on number of clusters used

Figure 11: Plot of within groups sum of squares depending on number of clusters used

Figure 12: Clusterplot of the 4 k-means clusters formed according to number of clusters promising least sum of squares and lowest complexity

Figure 12: Clusterplot of the 4 k-means clusters formed according to number of clusters promising least sum of squares and lowest complexity

Figure 13: Biplot of d18O and Sr/Ca values (all pike, all years) showing the life history strategy (anadromous in red, marine in blue and Freshwater in green) along with the assigned cluster (1-4, different shapes) based on k-means clustering.

Figure 13: Biplot of d18O and Sr/Ca values (all pike, all years) showing the life history strategy (anadromous in red, marine in blue and Freshwater in green) along with the assigned cluster (1-4, different shapes) based on k-means clustering.